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ABSTRACT 

An excess of sight-lines close to Lyman-break galaxies (LBGs) with little or no ab- 
sorption in QSO absorption spectra has been reported and has been interpreted as the 
effect of galactic winds on the Intergalactic Medium. We use here numerical simula- 
tions to investigate the flux probability function close to plausible sites of LBGs. We 
show that the flux distribution near our LBGs in the simulation depends strongly on 
redshift, and is very sensitive to the averaging procedure. We show that a model with- 
out galactic winds and a model with a wind bubble size of 0.5 h~^Mpc (comoving) are 
equally consistent with the new determination of the conditional flux distribution by 
Adelberger et al. (2005). Models with the larger bubble sizes ( >1 h~^Mpc) suggested 
by the previous observations of Adelberger et al. (2003) based on a much smaller sam- 
ple at higher redshift are not consistent with the new data. We, therefore, argue that 
the volume filling factor of galactic winds driven by LBGs may be much smaller than 
previously thought and that most of the metals responsible for the metal absorption 
associated with the low column density Lya forest are unlikely to have been ejected 
by LBGs. 
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1 INTRODUCTION 

Feedback from supernovae is an important ingredient in sce- 
narios of galaxy formation. Observations of high-redshift 
galaxies (Pettini et al. 1998; 2001), nearby starbursts (Heck- 
man et al. 1990) and dwarf galaxies (Heckman et al. 1995; 
2001) suggest that supernovae-driven winds are ubiquitous 
in star-forming galaxies. These could explain among oth- 
ers the high temperature of the high-redshift IGM (Cen & 
Bryan 2001), the faintness of the soft X-ray background 
(e.g.. Pen 1999) or the widespread detection of metals in 
the spectra of quasars (Cowie et al. 1995; Ellison et al. 1999; 
Schaye et al. 2000; Songaila 2001). 

In the standard picture, galactic winds are driven by 
the energy produced by supernova (SN) explosions of young 
massive stars. However, detailed modelling of a multi-phase 
ISM remains challenging (e.g. ©striker & McKee 1988; Ef- 
stathiou 2000) and predicting physical properties of galac- 
tic winds proves difficult. Simulations incorporating simple 
phenomenological prescriptions have succeeded in produc- 
ing SNe driven galactic winds (Mac Low & Ferrara 1999; 
Springel & Hernquist 2002) . However, the results of different 
numerical studies often do not agree (e.g.. Croft et al. 2002; 
Theuns et al. 2002), and the impact of winds on the IGM 



remains controversial. In particular, the volume fraction of 
the IGM affected by galactic winds, the typical mass of the 
galaxies responsible for the enrichment with the metals as- 
sociated with the low column density Lya forest, and the 
redshift at which the enrichment occurred remain a matter 
of debate (Gnedin 1998; Cen & Ostriker 1999; Aguirre et al. 
2001a, 2001b; Thacker, Scannapieco & Davis 2002; Bertone, 
Stoehr & White 2004; Haehnelt & Fieri 2004; Aguirre et al. 
2005; Fieri, Schaye & Aguirre 2005; Forciani & Madau 2005; 
Ranch et al. 2005). 

The Lyman-break technique (e.g. Steidel et al. 1996) 
has enabled the identification of large numbers of star- 
forming high-redshift galaxies. The large velocity shifts 
(100 — 1000 kms~^) between stellar and interstellar lines 
in the spectra of LBGs are strong evidence that these galax- 
ies drive powerful galactic winds (Fettini et al. 1998, 2001; 
Heckman 2000). However, these spectroscopic measurements 
do not tell us much about the location and volume filling fac- 
tor of the outflowing gas. Measurements of the transmitted 
flux in QSO absorption spectra with lines-of-sight close to 
galaxies can provide precious information on the distribution 
of ionized hydrogen in the surroundings of these galaxies. 
Some of this gas has presumably been shock-heated by out- 
flows. In a seminal study, Adelberger et al. (2003, hereafter 
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A03) measured the flux of the Lya forest in the absorption 
spectra of QSOs as function of the distance to the nearest 
LBG. At intermediate distances of 5-1 h~^Mpc (comoving) 
they found a decreeise of the mean flux level as expected 
if LBGs are located in overdense regions. The few (three) 
galaxies in their sample lying closer than 0.5 h~^Mpc to 
the line-of sight showed, however, a nearly complete lack 
of Lya absorption at the redshift of the galaxies. These re- 
sults sparked ofi^ a series of papers that attempted to ex- 
plain this by invoking feedback from ionizing radiation (e.g. 
Maselli et al. 2004) or from galactic winds (Croft et al. 2002; 
BruscoU et al. 2003; Kollmeier et al. 2003a, 2003b; Des- 
jacques et al. 2004; Kollmeier et al. 2005). AH these stud- 
ies had severe difficulties to reproduce the observed large 
decrease of the absorption at separation s ^0.5 h~^Mpc 
from the galaxies. Rather large bubble sizes (1 h~^Mpc or 
larger) and extremely efficient feedback due to SN (Croft 
et al. 2002) and/or high star formation rate (Maselli et al. 
2004) were required to explain the observations. It was also 
pointed out that thermal motions and coherent peculiar ve- 
locities strongly affect the galaxy-Lya absorption correlation 
at small separations (Kollmeier et al. 2003a). Desjacques et 
al. 2004 showed that these make it very difficult to explain 
the rather sudden increase of the mean flux level at a dis- 
tance of 0.5 h"^Mpc. 

Adelberger et al (2005,A05) have presented new results 
for the transverse proximity effect of LBGs in QSO absorp- 
tion spectra for a larger sample at somewhat lower redshift 
than A03. In the new sample, there are 24 galaxies at line- 
of-sight sepaxations smaller than 1 h~^Mpc. The sample is 
large enough to investigate the conditional probability dis- 
tribution of the flux at this distance. The picture is quite 
different to that of AOS where three out of three galaxies 
showed very little absorption. The absorption close to these 
24 galajcies covers the full range of flux levels but the PDF 
still shows a peak at low absorption. A05 have compared 
their measurements to numerical simulations of Kollmeier 
et al.. They found good agreement but an excess of objects 
with little absorption. They attributed this to the lack of 
galactic winds that were not modelled in the simulations of 
Kollmeier et al. We use here our own simulations with galac- 
tic winds similar to those presented in Desjaqucs et al. (2004, 
D04) to investigate the implications of the new observations 
for the size of the wind bubbles in our model. 



2 THE MODEL 

2.1 The Lya forest and the galeixy distribution 

We model the flux distribution of QSO absorption spectra 
as a function of the separation from the nearest LBG as 
described in D04. The modelling is based on a ACDM sim- 
ulation whose properties are discussed in detail in Stoehr et 
al. (2002). Note that we only consider the particle distribu- 
tion inside a box of size L = 30 h^^Mpc that belongs to a 
spherical region of high resolution. The particle mass in this 
high-resolution region is m = 1.66 x 10^ M©/h. We average 
the statistics of the simulated flux from three snapshots at 



redshift z = 3, 2.5 and 2, to facilitate a comparison with the 
observations, whose median redshift is z = 2.3. 

We associate LBGs with dark matter (DM) haloes iden- 
tified with a Friends-of-Friends group finding algorithm, as- 
suming that a DM halo contains not more than one LBG 
(e.g. Adelberger et al. 1998). Given that LBGs have a wide 
range of stellar masses, and that the relation between LBGs 
and DM haloes is still somewhat uncertain, we have inves- 
tigated two models as in D04. In the "staxburst" model, we 
assume that the observed LBGs reside in low mass haloes 
(e.g. Lowcnthal et al. 1997; Tragcr et al. 1997; Somerville 
et al. 2001). We randomly select DM haloes such that their 
number density in the simulation is similar to the comoving 
observed number density tIlbg of high-redshift LBGs. We 
take n^gQ to be constant in the redshift range 2 < z < 3, 
"lbg = 0-004 h^Mpc"-'* (Giavalisco & Dickinson 2001). In 
the "massive-halo" model, we assume that the high-redshift 
LBGs are the progenitors of the present-day massive and 
luminous galaxies (e.g. Steidel et al. 1996; Adelberger et al. 
1998; Bagla 1998; Haehnelt, Natarajan & Rees 1998). We 
pick a minimum mass Mmin = 5 x 10^^ M0/h such that the 
number density of haloes M > Mmin in the simulation is 
similar to the observed LBG density, tIlbq . 

Synthetic spectra are extracted from the DM simulation 
in the usual way. The spectra are normalized to reproduce 
the mean transmitted flux of observed high-resolution spec- 
tra, (F) = 0.68, 0.79 and 0.85 at 2 = 3, 2.5 and 2 respectively 
{e.g. McDonald et al. 2000). 



2.2 The effect of wind bubbles 

Our simple model for the effect of wind bubbles assumes 
that shocks produce long-lived fully ionised spherical bub- 
bles around galaxies (Croft et al. 2002; Kollmeier et al. 
2003a, 2003b; Weinberg et al. 2003; D04). We take the neu- 
tral hydrogen Hi density to be zero in a spherical region of 
radius R centered on the DM haloes which we have identified 
as hosts of LBGs. In the simulation, most of the haloes lie 
along filaments and the gas distribution is often anisotropic 
(e.g. Croft et al. 2002, Springcl & Hcrnquist 2003). Theuns 
et al. (2002) have also shown that winds propagate more 
easily into the low density IGM than in the filaments. We 
therefore caution that our assumption of a spherical wind 
bubble will only be a reasonable approximation for strong 
isotropic winds. Furthermore wind bubbles do not necessar- 
ily live for a Hubble time. 



3 COMPARISON WITH OBSERVATIONS 

3.1 Cosmic variance and galaxy redshift 
uncertainty 

The new sample of A05 is significantly larger than that of 
A03. At separation s < 1 h^^Mpc, there arc 24 galaxies with 
good redshift determination. Even though the new sample 
is larger than the AOS sample, it is still crucial to assess the 
statistical errors on the derived flux probability distribution 
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Figure 1. The mean flux F as a function of tlie distance s to 
the nearest halo in the starburst model for different radii of the 
wind bubble, 0.5, 1 and 1.5 h~^Mpc (top panel). The solid curve 
shows F for a model without galactic winds {R = 0) . The horizon- 
tal dashed line shows the mean transmissivity (/)=0.765 of the 
observed sample. The shaded area is our cosmic variance estimate 
for the model without galactic winds calculated from the cumu- 
lative distribution A^gal shown as the light (blue) histogram at 
the bottom. The filled symbols (top) and the dark bars (bottom) 
show F and A^qal of the observed sample of A05. 



resulting from the finite number of galaxies and QSO spec- 
tra. We focus on the conditional PDF P[D, s < 1) of the flux 
decrement D = 1 — F a.t a, separation s < 1 h~^Mpc, and 
the conditional transmitted flux F{s). We estimate the sta- 
tistical errors using Monte-Carlo realizations of the flux dis- 
tribution. We found that a sample of about 80 LOS through 
our simulation box gives a cumulative "galaxy" function 
Ngai^[< s) close to that of the observed sample at sepa- 
ration s < 1 h~^Mpc (see Fig. 1). The simulated sample has 
28 pairs with separation s < 1 h~^Mpc, while the observed 
sample has 24 pairs on that scale. The discrepancy at large 
scales is due to the (small) size of the observed fields. The 
number of haloes with mass M > 10" Mo/h in the simu- 
lation is about 1000. We should thus not have oversampled 
the halo catalogue too much. To account for redshift errors, 
we add a Gaussian error with variance An — 60 kms~^. 
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Figure 2. The differential probability distribution of the flux 
decrement B at a separation s < 1 h~^Mpc from a LBG, in the 
starburst and massive-halo models (top and bottom respectively) . 
P(D\s < 1) is plotted for the models of Fig. 1. The histogram 
shows the measurements of AOS. The vertical bars indicate the 
cosmic variance in bins of AD = 0.1 for the model R = 0. 



the cumulative halo distribution plotted as the fight (blue) 
histogram in the bottom panel. On scales s <A h~^Mpc our 
estimate of cosmic variance is consistent with that of A05, 
though somewhat larger by ~20 per cent. Fig. 1 shows that 
the 7? = and R — 0.5 models are consistent with the data 
of A05, while the models with R = 1 and 1.5 significantly 
overestimate F at separation s <1 h~^Mpc. 



3.3 The conditional probability distribution 

We calculate P{D\s < 1) and its cosmic variance error using 
Monte-Carlo realizations of 28 sight lines, where D is the 
mean fiux decrement of all pixels within 1 h~^Mpc from a 
LBG. We assume that there is one sight line per LBG. We 
normalize P{D\s < 1) such that the observed PDF gives the 
observed number of galaxies in each bin. 



3.3.1 The conditional PDF in the starburst model 



3.2 The conditional transmitted flux 

Fig. 1 compares the conditional fiux of our starburst model 
with the measurements of AOS, shown as filled symbols. The 
assumed bubble radius is i? = (solid), 0.5 (long dashed), 1 
(short dashed) and 1.5 h~^Mpc (dot), respectively. F has 
been computed from the full three-dimensional distribution 
of the fiux distribution in the simulation. Results obtained 
from the various snapshots have been rescaled to the mean 
transmissivity (/) = 0.765 of the observed sample (shown 
as the horizontal dashed line) before averaging. The shaded 
area shows the la error for the R = model computed from 



In Fig. 2 we compare the conditional PDF of the models 
plotted in Fig. 1 with the measurements of A05. The narrow 
vertical bars indicate the la error (in bins of AD — 0.1 for 
the model with 7? = 0). All models show a broad distribu- 
tion whith peaks at small and large flux decrements. The 
Kolmogorov-Smirnov (KS) probability that the observed 
distribution is drawn from the same PDF as our simulated 
distribution is highest for the model without galactic winds 
(0.315). It is still reasonably high for the model with R — 0.5 
(0.130), but decreases sharply for the models with larger 
bubble radii ( ;^10~''). The rather large errors suggest that 
statistical fluctuations could be responsible for the lack of 
LBGs with flux decrements in the range D ~ 0.3 — 0.6 as 
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well as the peak at D — 0.65, which is roughly 3cr above 
the mean. The existence of LBGs with D < 0.2 in the 
R = model may be surprising. One may have expected 
most galaxies to have D as they reside in dense re- 
gions with significant amounts of neutral hydrogen. Note, 
however, that the Hi density decreases rapidly away from 
LBGs. In fact, we found that ^55 per cent of the LBGs 
with D < 0.2 have a impact parameter > 0.8 h~^Mpc. This 
fraction is only ~ 35 per cent averaged over all flux decre- 
ments. Lines-of-sight with weak absorption have preferen- 
tially a large impact parameter in the R = model. Note 
that this trend appears not to be present in the observa- 
tions, where only 25 per cent of the LBGs with D <0.2 have 
an impact parameter > 0.8 h~^Mpc (Table 3 of A05). 



3.3.2 The conditional PDF in the massive-halo model 

The sensitivity of P{D\s < 1) to the bubble radius R sug- 
gests that it may also depend on the properties of the haloes 
hosting the LBGs. The bottom panel of Fig. 2 demonstrates 
that this is indeed the case. The conditional PDF is shown 
for our massive-halo model with the same bubble radii as for 
the starburst model. Again, all models have a broad distribu- 
tion with peaks at small and large flux decrements similar to 
that of the starburst model. However, the fraction of LBGs 
with flux decrement D < 0.2 is somewhat smaller. This is 
most likely due to the increase of the typical infall veloc- 
ities with increasing halo mass (D04). In the massive-halo 
scenario, the KS probability is highest for the model with 
bubble radius R = 0.5 h"^Mpc (0.151). The probability is 
only 0.014 in the R — model, and <^ 0.01 for the models 
with R> 1 h~^Mpc. A significant fraction of ionised bubbles 
with R ~ 0.5 h~^Mpc is thus required in the massive-halo 
scenario to reproduce the observed distribution. Note, how- 
ever, that the rather large errors make the interpretation 
of the KS probabilities problematic. Our results are con- 
sistent with our previous study where we also found that 
models where LBGs are mainly starbursts in small mass 
haloes matches the observations with smaller bubble radii 
than models where massive haloes host the LBGs (D04). 



3.3.3 Sensitivity to redshift and averaging procedure 

The conditional probability of the models shown in Fig. 2 
is an average of three equally weighted snapshots at z = 3, 
2.5 and 2. However, P{D, s < 1) evolves noticeably over that 
redshift range. This is shown in the top panel of Fig. 3, where 
P{D, s < 1) is plotted as a function of redshift for the star- 
burst model with _R = 0. Measurements at z = 2 contribute 
most to the low flux decrement tail of the PDF. The pre- 
dicted mean number of galaxies with D < 0.2 is n<o.2 ~ 2.6, 
5.0 and 8.6 in a sample with 24 objects sd, z — 3, 2.5 and 
2, respectively. In the massive-halo model, these values drop 
to n<o.2 = 1.7, 3.4 and 5.1, respectively. The strong depen- 
dence of P{D, s < 1) on redshift reflects the evolution of the 
unconditional PDF of the Lya flux in that redshift interval. 
The selection function of the A05 sample peaks at z — 2.3. 
Among the 24 galaxies with impact parameter < 1 h~^Mpc, 
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Figure 3. The differential probability distribution at separation 
s < 1 h^'^Mpc from a LBG as a function of the flux decrement 
D in the starburst model without galactic winds (i? = 0). The 
top panel shows the redshift dependence, and the bottom panel 
demonstrates the sensitivity to the averaging procedure (see text 
for details). The histogram is the data of AOS. 



only 5 have a redshift z > 2.5. We have also tried a weight- 
ing which matches the observed redshift distribution more 
closely and got similar results. 

The bottom panel of Fig. 3 demonstrates that the shape 
of the conditional PDF is also strongly sensitive to the av- 
eraging procedure used to calculate D. The solid curve is 
for our default assumption, and is computed as in A05. To 
obtain the long dashed curve, we averaged over four lines- 
of-sight per LBG. The short dashed curve shows the PDF 
obtained if pixels in the spectrum are considered individu- 
ally. Averaging over four lines-of-sight per LBG moves the 
flux decrement in the tails of the distribution closer to the 
mean. When all pixels are considered separately the peaks at 
small and large flux decrement are enhanced. These findings 
are consistent with the results of KoUmeier et al. (2003a). 



4 SUMMARY AND DISCUSSION 

We have revisited the proximity effect of LBGs on the ob- 
served flux distribution in QSO absorption spectra and have 
investigated the implications of the new data of A05 for the 
bubble size and volume filling factor of galactic winds of 
LBGs. We have used mock spectra calculated from Lya for- 
est simulations with the galactic wind prescription of D04 to 
model the conditional probability distribution of the flux as 
a function of the line-of-sight distance to the closest LBG. 
As in D04 we have considered two scenarios for the LBGs, 
a starburst model and a massive-halo model. We have used 
Monte-Carlo realizations of the flux distribution to assess 
the statistical uncertainties. Our main results are as follows. 
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• Our simulations give a broad probability distribution 
of the flux for lines-of-sight passing close to LBGs (s < 
1 h~^Mpc ) with a peak at small and large flux decrements 
consistent with that observed by A05 within the expected 
errors for a sample of this size. 

• None of the models leads to the observed lack of gala^x- 
ies with decrement 0.3 < D < 0.6, or the observed excess 
at 5 ~ 0.7. However, error bars are laxge and statistical 
fluctuations may be responsible for these features. 

• For the starburst model the conditional PDF of the 
model with no galactic winds agrees best with the observa- 
tions, but the model with a bubble radius R = 0.5 h~^Mpc 
is also consistent. For the massive halo model the model with 
R = 0.5 h^^Mpc agrees best with the observation. 

• The conditional PDF of the flux depends strongly on 
redshift and the averaging procedure applied. 

Our results appear to disagree somewhat with the pro- 
diction of the SPH simulation without galactic winds shown 
in Fig. 15 of the A05 paper. Note, however, that the latter 
was obtained from a single output at z = S rescaled to the 
mean redshift of the observed sample. The strong depen- 
dence of P{D, s < 1) on redshift could explain the absence 
of the bimodality in the PDF of the SPH simulation. Of 
course, our simplified treatment of the IGM physics may 
also be the culprit for the difference. 

The good agreement of our staxburst model without 
galactic winds with the meeisurements of A05 suggests that 
the volume fraction of the IGM affected by winds may be 
very small. A model where LBGs drive outflows extending 
substantially less than 0.5 h~^Mpc (comoving) can still ex- 
plain the peak in the conditional PDF at small flux decre- 
ment, and reproduce the mean flux level as a function of the 
distance to the closest LBG, irrespective of the exact prop- 
erties of the DM haloes hosting LBGs. Note that in such a 
model galactic winds from LBGs, while still being respon- 
sible for the observed strong CIV absorption close to these 
galaxies (A05), would nevertheless have very little effect on 
global statistical properties of the Lya forest (e.g. Theuns 
et al. 2002; Desjacques et al. 2004; McDonald et al. 2005). 
They would also not be responsible for the majority of the 
weaker CIV absorption which has a much larger volume fill- 
ing factor. 
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